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Trilayer graphene in the fractional Quantum Hall Effect regime displays a set of unique interaction- 
induced transitions that can be tuned entirely by the applied bias voltage. These transitions occur 
near the anti-crossing points of two Landau levels. In a large magnetic field (> 8 T) the electron- 
electron interactions close the anti-crossing gap, resulting in some unusual transitions between differ- 
ent Landau levels. For the filling factor u — ^, these transitions are accompanied by a change of spin 
polarization of the ground state. For a small Zeeman energy, this provides an unique opportunity 
to control the spin polarization of the ground state by fine tuning the bias voltage. 



Dirac fermions in monolayer and bilayer graphene with 
their remarkable electronic properties have received ex- 
traordinary scrutiny in recent years 0, Si- In an ex- 
ternal magnetic field, these systems exhibit unconven- 
tional quantum Hall effects [j, Q that are direct mani- 
festations of their rather unusual band structures 0j ■ 
As a consequence, the Landau level (LL) energies of 
these systems are very different from those of conven- 
tional two-dimensional electrons systems (2DESs). More 
specifically, in monolayer graphene the Landau level en- 
ergies exhibit a square root dependence on the applied 
field Q, while in bilayer graphene one finds a linear de- 
pendence 0. On the other hand, interactions among 
Dirac fermions in the fractional quantum Hall effect 
(FQHE) n^.[Tl| regime reveals several rather unexpected 
and intriguing effects in monolayer [13] and bilayer [13] 
graphene. Recent experimental observation of the FQHE 
in monolayer graphene [3, [13 have indeed confirmed 
the important role electron-electron interactions play in 
these systems. Clearly, the dynamics of Dirac fermions 
are sensitive to the number of graphene layers present 
in the system and their stacking arrangements. Quite 
expectedly, the attention has now shifted to the investi- 
gation of the electronic properties of Dirac fermions in 
trilayer graphene (TLG). 

A TLG consisting of three coupled graphene layers has 
a very unique electronic energy spectrum. Within the 
nearest-neighbor inter-layer coupling approximation the 
energy spectrum of TLG with Bernal stacking consists 
effectively of decoupled single-layer graphene and the bi- 
layer graphene energy spectra. Therefore the TLG al- 
lows us to study the energy spectra of both the massless 
and massive Dirac fermions within a single system. In a 
strong perpendicular magnetic field, the LL energy spec- 
trum of TLG becomes a combination of Landau levels of 
single-layer and bilayer graphene [l6| . The spectrum ex- 
hibits many crossings of the Landau levels as a function 
of the magnetic field. At the crossing points the Landau 
levels are highly degenerate. The degeneracy is lifted 
when higher-order inter-layer coupling terms are taken 




FIG. 1: (a) Schematic illustration of the ABA stacking of 
TLG. Each graphene layer consists of two inequivalent sites 
A and B. The inter-layer and intra-layer hopping integrals, 7^, 
show the couplings, which are included in the single-particle 
Hamiltonian ([1]). (b) Landau level energy spectrum of TLG 
in a magnetic field of 15 T as a function of the bias voltage, U. 
Two red lines show anti-crossing for U ~ 400 meV. The corre- 
sponding LLs are labelled as level-1 and level-2, respectively. 
The LLs 1 and 2 belong to the set of Landau levels with pa- 
rameter n — 0. The inset shows the region of anti-crossing. 
The anti-crossing gap is « 2.6 meV~30 K. 



into account, resulting in several unusual properties of 
the quantum Hall effect in trilayer graphene [17[. 

While the single-particle features of the quantum Hall 
effect are interesting, here we show that by introducing 
interactions among the Dirac fermions in a TLG in the 
FQHE regime, we witness several very unique properties 
of the TLG that goes far beyond the mere level cross- 
ings observed in the integer QHE. We found several LL 
repulsions and level crossings which resulted in some in- 
teresting spin transitions among the LLs in this system 
that have no analogues in the interaction-induced spin- 
reversed ground states and elementary excitations dis- 
covered earlier in conventional electron systems [l^, [l^ . 
These spin transitions in the TLG are driven by an ap- 
plied perpendicular bias field for a fixed magnetic field, 
and therefore we expect these novel transitions to be en- 
tirely tunable. 

In what follows, we only consider the Bernal or ABA- 
stacking of our TLG. In the tight-binding approximation 
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the Hamiltonian of TLG is characterized by the intra- 
layer hopping integral, Jq = 3.1 eV, and inter- layer hop- 
ping integrals, 7^ = 0.39 eV, 72 = -0.028 eV, 74 = 0.041 
eV, and 75 = 0.05 eV, corresponding to different types 
of the inter-layer coupling, shown schematically in Fig. 
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netic field, the Hamiltonian of a TLG for a single valley, 
e.g., valley K, takes the form (TgI. [iTj 
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where Vq = (\/3/2)a7Q/?i sa 10^ m/s, '(;4 = {^/3/2)a'-f^/h, and ttj. = tt^ ± tt^^. Here n = p + eA/c is the generalized 
momentum. The parameter S = 0.046 eV is the difference between the on-site energies of two sublattices within a 
single graphene layer [l3|- The bias voltage, U, is introduced in the Hamiltonian ([T]) as the potential difference, i.e., 
the on-site energy difference, between layers 1 and 3. Here we assume that the potential of layer 1 is zero, while the 
potential of layer 2 and 3 are U/2 and U, respectively. The bias voltage, U, is considered as the parameter of the 
system, which can be varied externally. 



The LLs of a TLG can be obtained from the Hamil- 
tonian matrix ([1]). The corresponding wave functions 
are parametrized by the integer n and can be expressed 
through the conventional (non-relativistic) Landau level 
wave functions, 0„ ^ as 

C'l0n-|-2,m ^ 
C'20n-|-l,m 

, (2) 

where m is the intra-Landau level parameter, e.g., the 
angular momentum, and Cj are constants. Therefore, 
the LL wave functions of a TLG are combinations of n, 
n + 1, and n + 2 non-relativistic Landau functions. 

In what follows, we consider a many-electron system 
partially occupying a single LL of the TLG. We study the 
properties of these systems in the FQHE regime, specif- 
ically for the filling factors, v = |, |, and |. In these 
cases, the conventional non-relativistic system shows in- 
compressible behavior with a finite energy gap [l^, [HI ■ 
The interaction properties of the many-electron system 
occupying a single Landau level are completely deter- 
mined by the Haldane pseudopotentials Vm'^ [l^l, which 
are the interaction energies of two electrons with relative 
angular momentum, m. They are determined by mak- 
ing use of the LL wave functions ^ '^ll. We numeri- 
cally evaluate the FQHE state in a TLG by considering a 
finite-size system of N electrons in a spherical geometry 
[20j with interaction potentials determined by the Hal- 
dane pseudopotentials. The radius of the sphere is \/S£q, 
where 25' is the number of magnetic fluxes through the 



sphere in units of the flux quantum, and £q = -y/ eh/cB is 
the magnetic length. The parameter S also determines 
the number of single-particles states, 25-1-1, and for flnite 
number of electrons - the flUing factor of the system. 

From the Hamiltonian ([Ij we numerically evaluate the 
single-particle LL energy spectrum. The TLG LLs are 
parametrized by the integer n [see Eq. For each 

n there are 6 LLs in a TLG. A typical LL spectrum 
is shown in Fig. 1 (b). The spectrum as a function of 
the bias voltage (or the magnetic held) shows crossing 
and anti-crossing of the energy levels. The anti-crossing 
gap, shown in Fig. 1 (b), is around 2.6 meVw 30 K for 
a magnetic field of 15 Tesla. Below we show that the 
FQHE in TLG has non-trivial and unique interaction- 
induced properties near these anti-crossing points. We 
study the behavior of the system near the special anti- 
crossing point shown as inset in Fig. 1 (b). This point 
corresponds to anti-crossing of the TLG LLs with n = 0. 
We label the corresponding levels as Landau level- 1 (LL- 
1) and Landau level-2 (LL-2) [Fig. 1 (b)] and consider the 
FQHE states only in these levels. For each FQHE state 
in the corresponding LL-1 and LL-2 we have evaluated 
the ground state energy per particle and the excitation 
gap for the incompressible states. 

We first consider the fundamental v = i-FQHE states 
in LLs 1 and 2, respectively [Fig. 2]. The system be- 
haves very differently for weak and strong magnetic fields. 
In a weak magnetic field (i? = 5 T) the many-particle 
states show anti-crossing (Fig. 2a) similar to the single- 
particle levels [Fig. 1 (b)]. This anti-crossing is clearly 
visible in the dependence of the i-FQHE gaps since the 
values of the FQHE gaps in LL-1 and LL-2 are inter- 
changed when the system goes through the anti-crossing 
point. The system shows an interesting behavior exactly 
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FIG. 2: Numerical results for finite-size v = |-FQHE systems 
in LL 1 and LL 2, shown as a function of the bias voltage near 
the anti-crossing points. The magnetic field is 5 Tesla [panels 

(a) and (c)] and f5 Tesla [panels (b) and (b)]. Panels (a) and 

(b) show the ground state energy per particle, while panels 

(c) and (d) show the excitation gaps at the corresponding 
energy levels with filling factor i/ — ^. Blue arrows in panel 

(d) indicate jumps of the FQHE gaps at the level crossing 
in panel (b). The number of electrons is A'' = 9 and the 
parameter of the sphere is 25 = 24. 

at the anti-crossing point. Here, due to a mixture of 
the single-particle wavefunctions of LL-1 and LL-2, the 
many-particle interaction properties are enhanced in LL- 
1 while suppressed in LL-2. As a result, at the anti- 
crossing point, the -l-FQHE gap in LL-1 has a maximum 
while the i-state in LL-2 becomes compressible with a 
vanishing gap [Fig. 2 (c)]. Due to a larger cohesive energy 
of the incompressible state compared to the compressible 
one, the many-particle anti-crossing gap shows a small 
enhancement relative to the single-particle value by 0.4 
meV w 5 K. Experimentally, the anti-crossing proper- 
ties of TLG in a small magnetic field can be observed 
by studying the FQHE in LL-2. In such a system, with 
increasing bias voltage one would observe a transition 
FQHE - no FQHE - FQHE within a single LL, just as 
we predicted earlier for bilayer graphene |13]. 

Near the anti-crossing point in a large magnetic field 
{B — 15 T), TLG shows several novel features [Fig. 2 
(b,d)]: The anti-crossing of the single-particle energy 
levels becomes double crossings for the many-particles 
states. This means that the cohesive energy of the many- 
particle state in LL-2 is larger than that in LL-1, and 
this difference overcomes the anti-crossing gap. The rea- 
son for such a behavior is the change in the interac- 
tion strength in LL-1 and LL-2. For a large magnetic 
field the many-particle interaction potential at the anti- 
crossing point becomes stronger in LL-2 and weaker in 
LL-1, which is opposite to what we see for a weak mag- 
netic field [Fig. 2 (a,c)]. As a result, the FQHE gap in 
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FIG. 3: The anti-crossing gap corresponding to the anti- 
crossing of LL-1 and LL-2 [Fig. 1 (b)], versus the magnetic 
field. For different magnetic fields, anti-crossings occur for dif- 
ferent bias voltages. The black line and the squares describe 
a single-particle system, while the red line and the circles cor- 
respond to the many-particle f — i-FQHE system (Fig. 2). 
The single-particle anti-crossing gap is closed for the many- 
particle system for B ~ 8 T. The number of electrons in the 
many-electron system is N — 9 and the parameter of the 
sphere is 2S = 24. 



LL-2 has a maximum at the anti-crossing point, while the 
gap in LL-1 is suppressed. Therefore at the anti-crossing 
point the — -l-many-particle system is incompressible 
in LL-2 and compressible in LL-1. Since the incompress- 
ible FQHE state has a lower binding energy than that of 
the compressible state, this energy difference is enough 
to close the anti-crossing gap. Therefore, in a large mag- 
netic field and as a function of the bias voltage, we should 
expect the following behavior: The FQHE system, which 
initially for a small bias voltage, U < 400 meV, is in 
LL-1, occupies LL-2 at the anti-crossing point, U — 400 
meV, leaving LL-1 empty. With further increase of the 
bias voltage, U > 400 meV, the system returns to LL-1, 
while LL-2 becomes empty. These transitions between 
different LLs at the anti-crossing point are accompanied 
by jumps in the value of the FQHE gap as illustrated by 
blue lines in Fig. 2 (d). 

The strength of the FQHE, i.e., the magnitude of 
the FQHE gap, and correspondingly the cohesive energy 
of the FQHE states, is determined by the short-range 
properties of the interaction potential, i.e, the Haldane 
pseudopotentials at small values of the relative angu- 
lar momentum, m. Therefore, at a weak magnetic field 
the short-range interaction strength at the anti-crossing 
point is enhanced in LL-1 (lower energy level), while 
for high magnetic fields the interaction strength at the 
anti-crossing point is increased in LL-2. This results in 
a weak enhancement of the many-particle anti-crossing 
gap for weak magnetic fields and strong suppression of 
the many-particle anti-crossing gap for strong magnetic 
fields. In Fig. 3, we present the anti-crossing gaps for 
single-particle and many-particle — i-FQHE systems 
for different magnetic fields. For small values of the 
magnetic fields, 1.5T < B < 8T, the many-particle gap 
clearly shows an enhancement compared to that for the 
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single-particle case. For _B 8 T the many-particle anti- 
crossing gap closes, and for i? > 8 T the anti-crossing in 
a single-particle system becomes a double-crossing in the 
many-particle FQHE system. 

It is known that the FQHE ground state of conven- 
tional semiconductor systems can be spin-polarized or 
spin-unpolarized. The polarization properties of the sys- 
tem are determined by the filling factor and strength of 
the applied magnetic field [ll|, 13| • While in conven- 
tional systems the = ^ FQHE state is always fully 
spin-polarized, for filling factors v = ^ and v — there 
is a competition between the energies of spin-polarized 
and spin-unpolarized incompressible states [18|. With 
increasing strength of the magnetic field, the Zeeman en- 
ergy of electrons favors the spin-polarized state, which re- 
sults in possible spin-transitions in the system by varying 
the magnetic field. Those theoretical predictions subse- 
quently received experimental confirmation [19,. i22ri25| . 
The TLG system shows a different type of spin tran- 
sitions realized at the anti-crossing points, which could 
also be probed experimentally. 

We have analyzed the spin properties of the FQHE 
states in TLG for filling factors v = ^ and = |. 
In Fig. 4 the results for v = |-FQHE state are shown 
for LL-1 and LL-2 without including the Zeeman en- 
ergy. The black and red lines in Figs. 4 (a,b) correspond 
to spin-polarized and spin-unpolarized systems, respec- 
tively. The general behavior of the system is similar to 
that of the v = i-FQHE state. In a small magnetic field 
the system shows stronger interactions in LL-1, while for 
larger magnetic fields the anti-crossing of energy levels 
becomes double crossings. For a small magnetic field, 
the ground state of the v — |-FQHE system is mainly 
spin-polarized with only a small region of bias voltages, 
U , when the system becomes spin-unpolarized in LL-2. 
Therefore, for a weak magnetic field, B ^ 5 T, the = |- 
FQHE system in LL-2 should show spin transition into 
an unpolarized state within a narrow interval of U at the 
anti-crossing point. A strong Zeeman energy will how- 
ever suppress this spin transition. 

In a large magnetic field [Fig. 4 (b)], the v = |-FQHE 
system shows interesting spin properties. While in LL- 
1 the y = |-FQHE ground state is spin-unpolarized for 
small values of J7, it becomes spin-polarized at a large 
bias voltage, U > 420 meV. In LL-2 the ground state 
is spin-unpolarized for all values of U. Finally, combin- 
ing these two behaviors and comparing the ground state 
energies of different systems [Fig. 4 (b)], we predict the 
following novel spin transitions. If the system is initially 
in LL-1 for the = |-FQHE state, then with increas- 
ing bias voltage the system will undergo the following 
transitions: spin-unpolarized state in LL-1 spin- 
polarized state in LL-2 spin-unpolarized state in 
LL-1 <^ spin-polarized state in LL-1. What is remark- 
able here is that, spin polarization of the ly = |-FQHE 
system in TLG can be controlled by fine tuning the bias 
voltage ~ a possibility that never existed in the FQHE 
regime of conventional systems. We have also studied 
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FIG. 4: The ground state energy per particle [panels (a) and 
(b)] and the excitation gaps [panels (c) and (d)] for v — ^- 
FQHE in LL-1 and LL-2 versus the bias voltage near the 
anti-crossing points. The excitation gaps are shown only for 
the spin-polarized systems. The magnetic field is 5 Tesla 
[panels (a) and (c)] and 15 Tesla [panels (b) and (b)]. The 
black and red lines in panels (a) and (b) correspond to spin- 
polarized and spin-unpolarized systems, respectively. For the 
spin-polarized system, the number of electrons is A'^ = 16 and 
the parameter 25' is 24, while for the spin-unpolarized system 
N = 10 and 25" = 12. 



the 1/ = I filling factor in LL-1 and LL-2. The gen- 
eral properties of the anti-crossing is similar to those of 
the v — i-FQHE system, i.e., an enhancement of the 
anti-crossing gap for a small magnetic field and strong 
suppression of the anti-crossing gap for a large magnetic 
field. However, for both LL-1 and LL-2 the z/ = |-FQHE 
ground state is spin-polarized, which excludes the possi- 
bility of spin transitions in that system. 

In conclusion, trilayer graphene exhibits several unique 
electronic properties near the anti-crossing points of two 
Landau levels. In the FQHE regime, the electron-electron 
interaction strongly renormalizes the anti-crossing gap. 
In a weak magnetic field (B < 8 T), the many-body in- 
teraction enhances the anti-crossing gap, resulting in a 
non-monotonic dependence of the excitation gaps on the 
bias voltage. In a large magnetic field (B > 8 T), the 
electron-electron interaction strongly suppresses and fi- 
nally closes the anti-crossing gap. In such large magnetic 
fields, the spin-polarized FQHE system shows nontrivial 
transitions as a function of the bias voltage, which are ac- 
companied by jumps of the FQHE excitations gaps. In a 
large magnetic field {B > 8 T) the TLG displays unique 
spin polarizations with controllable spin transitions. By 
varying the bias voltage, the v = |-FQHE system can be 
switched from spin-polarized to spin-unpolarized states. 
Various experimental techniques have been developed in 
the past to study spin transitions in the FQHE regime for 
conventional electron systems. These include measure- 
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ments in a tilted magnetic field [l9[ , application of hydro- 
static pressure [12], optically pumped nuclear magnetic 
23|, resistively detected nuclear magnetic res- 
time- resolved radiative recombination 12511 . 



resonance 
onance [Ij] 

among others. Similar experimental studies in trilayer 



graphene will undoubtedly uncover a wealth of informa- 
tion about charge and spin transitions revealed by a TLG 
in the FQHE regime. 
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